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Abstract —A design of rate-compatible polar codes suitable for 
HARQ communications is proposed in this paper. An important 
feature of the proposed design is that the puncturing order is 
chosen with low complexity on a base code of short length, which 
is then further polarized to the desired length. A practical rate¬ 
matching system that has the flexibility to choose any desired 
rate through puncturing or repetition while preserving the 
polarization is suggested. The proposed rate-matching system is 
combined with channel interleaving and a bit-mapping procedure 
that preserves the polarization of the rate-compatible polar code 
family over bit-interleaved coded modulation systems. Simula¬ 
tion results on AWGN and fast fading channels with different 
modulation orders show the robustness of the proposed rate- 
compatible polar code in both Chase combining and incremental 
redundancy HARQ communications. 

L Introduction 

Polar codes, introduced by Arikan, are the first class of 
error-correcting codes that provably achieve the symmetric ca¬ 
pacity of memoryless channels with low-complexity encoders 
and decoders m. Arikan applied a linear transformation over 
GE(2) on an input vector of length N = 2^ bits based on the 
n-th Kronecker power F®" of the 2x2 matrix 


which is referred to as a kernel matrix. As the number of 
polarization levels n increases, the input channels are trans¬ 
formed into bit-channels that tend to become either noiseless 
or completely noisy channels under successive cancellation 
(SC) decoding. The fraction of the noiseless bit-channels 
approaches the symmetric capacity of the transmission channel 
m. The encoding and SC decoding complexities of Arikan’s 
polar codes are relatively low at N log 2 N. 

To design hybrid automatic repeat request (HARQ) rate- 
compatible (RC) codes for wireless systems, it is required 
that the information block length is fixed across all codes 
with different rates. It is also required that the transmitted 
bits can be flexibly chosen according to the desired code rate, 
and HARQ scheme. Chase combining (CC) or incremental 
redundancy (IR). Although there have been attempts to con¬ 
struct codes with different rates or lengths, they do not satisfy 
the flexible HARQ design requirements. The block lengths of 
Arikan’s polar codes were restricted to a power of two. Polar 
codes with arbitrary length can be constructed by replacing the 
kernel matrix F in ([T]) with a non-singular binary matrix of 
size f X £ ||2|- However, adjusting the polar code block lengths 
through changing the size of the kernel matrix implies different 


pairs of encoders and decoders for each block length, as well 
as increased decoder complexity at larger ^. Puncturing polar 
codes to obtain length-compatible codes has been viewed as a 
process of reduction of the size of the matrix F®", where the 
reduced matrix with the largest polarizing exponent is found by 
exhaustive search for each desired rate 13. Changing the code 
rate through information shortening and random puncturing of 
the output bits were investigated Another scheme changes 
the transmission rate through quasi-uniform puncturing as 
well transmissions of the input information bits 0. Several 
criteria based on minimizing the error probability have been 
investigated to chose a puncturing pattern 0 , where the search 
criterion is tested for each possible puncturing pattern, at each 
desired code length and desired code rate. 

One disadvantage of the previous algorithms is that they 
do not necessarily result in a HARQ RC family of codes. 
Moreover, the puncturing pattern has to be exhaustively found 
at each code length and code rate. In this paper, we propose 
an efficient algorithm that finds the best order for puncturing 
the output bits on a base polar code, while guaranteeing 
that the codes with different rates are nested. Our simulation 
results show that this scheme achieves a decoding performance 
similar to that when the optimal puncturing pattern is selected 
for a given rate and block length. We show how a regular 
puncturing pattern of a longer polar code can be derived 
from the puncturing pattern of the base polar code, while 
preserving the code polarization, based on the compound polar 
code construction 0. This reduces the complexity of the 
code design, where we do not need to exhaustively search 
for the puncturing at each desired code rate and code length. 
Moreover, the time and space complexities of operation at 
both the transmitter and the receiver are reduced. We propose a 
simple rate-matching scheme which implements the compound 
punctured polar codes and allows for flexible puncturing or 
repetition of code bits. Our scheme supports both HARQ 
Chase combining, when the bits of the retransmissions are the 
same as those of the first transmission, and HARQ incremental 
redundancy, when the retransmissions constitute some coded 
(redundancy) bits that have not been transmitted before and 
may also constitute some of the previously transmitted bits. We 
show how the proposed rate-matching structure also provides 
channel interleaving and bit-mapping for bit-interleaved coded 
modulation (BICM) while maximizing code polarization. 

The rest of the paper is organized as follows. Preliminaries 
of polar codes and estimation of their bit-channel error prob¬ 
ability are reviewed in Section |II] In Section [Till we describe 




the different aspects for the design of the proposed rate- 
compatible polar codes. Analysis of the achievable rates and 
the decoding performances of the proposed rate-compatible 
polar codes with different HARQ schemes are presented in 
Section |IV] Conclusions are made in Section FV] 

II. Preliminaries 

The polar code construction of length = 2” is based 
on the n-th Kronecker power F®" of the 2x2 matrix in 
O- Assume the information sequence is denoted as = 
{ui,U 2 , • ■ •, un) and the corresponding codeword is denoted 
as Xi = {xi,X 2 , ■ ■ ■ ,X]sf), where Xi, ut G {0,1} for 1 < 
i < N. The encoding of a polar code can be described by 
x^ = Ui where Bjsi is a bit-reversal permutation 

matrix m. Xi is transmitted through N independent copies 
of a binary input discrete memoryless channel. The received 
sequence is denoted as = (yi, 1 / 2 , • ■ •, t/Af)- The combined 
channel, Wm, between and is described through the 
transitional probabilities 

N 

WN{y?\u^) = ) = n W{y,\x,). ( 2 ) 

i=l 

Through recursive channel splitting and channel combining 
operations, a bit-channel for each input information bit Ui is 
defined as m 

= P{yi 

= H i^^WN{yi\u^)- 

Based on these bit-channel models, Arikan proposed a recur¬ 
sive SC decoding algorithm m. 

For uniform channel input, the error probability (EP) of 
the polarized bit-channels under genie-aided SC decoding on 
erasure channels can be recursively calculated through their 
Bhattacharrya parameters fTl. For other channels, the bit- 
channel error probabilities can be estimated through Monte- 
Carlo simulations of the genie-aided SC decoder. In case of 
SC decoding on AWGN channels, the bit-channel EPs can be 
approximated with reasonable accuracy by density evolution 
and approximating the outputs at each SC decoding step with 
Gaussian random variables [Si. Assume the all-zero codeword 
is transmitted and the variance of the AWGN channel is cr^. 
We denote the logarithmic likelihood ratio of the received yi 
as L{yi), which is assumed to have probability distribution 
From the basic encoding and decoding structure 
of polar codes in FiglT] we can find E(L(ui)) and E(L(u 2 )) 
from L{yi) and L{y 2 ), where E denotes the expectation. The 
error probability of bit-channel Ui is estimated by HI 

PiS,)=QiVnL{^^])). ( 4 ) 

With different L{yi) and T(j/ 2 ), the Gaussian approximation 
(GA) can be written as 

E[Liu,)] = </>-'(! - (1 - 0(E(L(t/i))]))(l - <^(E(L(y 2 ))]))) 

(5) 
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Fig. 1. Basic encoding and decoding structure of Polar codes 


where 
cj){x) = 


E[L{u 2 )]=E[L{y^)]+E[L{y 2 )], (6) 

1 ifa;>0, 

1, if a; = 0. 


Through recursively applying Q and (|6ll on the basic 
encoding and decoding structure of a polar code, E{L{ui)) 
of each information bit Ui can be calculated from the received 
L{yi) for 1 < i < W. In case an output channel is punctured, 
then its variance is set to infinity. 


III. Proposed rate-compatible polar codes 

In this section, we describe the different aspects of the 
proposed HARQ rate-compatible polar code design. 

Assuming an {N, k) mother polar code of rate k/N, where 
k is the cardinality of the information set I. X can be chosen 
to constitute of the k bit-channels with the lowest estimated 
bit-channel error probabilities to minimize the union bound on 
the decoder block error rate (BEER) m, 

P{£)<J2PiS.). (7) 

iei 

Our goal is to design a family of nested RC codes, with 
rates k/N or higher as in {k/N, k/{N — 1),..., k/{N — m)} 
by puncturing coded bits, as well as below k/N by repeating 
some of the coded bits at the same transmission. Also, the 
transmission rate can be allowed to go above 1 by transmitting 
less than k bits if N—m < k. In this design, the polar encoding 
is not constrained to be systematic. 


A. Progressive puncturing algorithm 

Let the desired punctured code rate be k/{N — m). The 
exhaustive search approach ||6| estimates the information set 
EP (|7]) after puncturing m bits, for all (^) possible puncturing 
patterns, and chooses the one that minimizes (|7]l. Alternatively, 
for different applications, the design criterion d?) can be 
changed. Whereas exhaustive search will find the optimal 
puncturing pattern for that rate, it does not guarantee that the 
puncturing patterns of different rates are nested. 

To address this, we propose an efficient greedy algorithm 
to find the best puncturing patterns for RC polar codes. The 
Progressive Puncturing Algorithm (PPA) finds the puncturing 
pattern that punctures m + 1 bits such that this puncturing 
pattern minimizes the design criteria (|7]) and is constrained to 
contain the m previously punctured bits. Hence, it results in a 
nested family of codes, where the code of rate k/ {N — m + 1) 
is obtained from the code of rate k/{N — m) by transmitting 
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Fig. 2. Progressive versus exhaustive search puncturing of polar codes 


the mth punctured bit. The set of indices of the punctured bits 
is denoted by V which is also referred to as the puncturing 
pattern. Given a certain puncturing pattern V of size m, the 
(m + l)th punctured bit is selected from the remaining N — m 
non-punctured coded bits by testing the design criterion N—m 
times. Hence, to find the puncturing order for 1 < m < N, 
the design criterion needs to be tested N{N + l)/2 for the 
PPA, which is significantly less than the 2^ searches required 
with exhaustive search. 

In Fig. 121 the decoder BLER is estimated using Q, where 
the bit-channel EP is found by the GA on an AWGN channel 
after code puncturing. Assuming an (32,16) polar code de¬ 
signed at SNR = 3 dB, the performance of the punctured polar 
code by PPA is shown to overlap with that of the punctured 
polar code found by exhaustive search, for toG{0,4,6,10}. 

B. Two-step polarization and regular puncturing 

Although, the PPA has significantly smaller complexity than 
exhaustive search, it is not practical to repeat it for every code 
length and rate. Moreover, even the PPA algorithm becomes 
infeasible at very large block lengths as it requires estimating 
the bit-channel error probabilities. In our proposed scheme, the 
PPA is only run on a base polar code of short length, e.g. N' = 
32. The base polar code is then further polarized into the polar 
code with the desired length N. This two-step polarization 
approach results in a regular puncturing pattern on the longer 
polar code of length N, which is derived from the puncturing 
sequence of the base code of length N'. Another advantage 
of the two-step polarization approach is that it avoids storing 
puncturing patterns for each code length and code rate. 

Consider an {N, k) mother polar code of length V = 2" and 
a base polar code of length N' = 2^. The encoding structure 
with generator matrix 32 ^ (F®”) can be decomposed into two 
stages as 82 ^ The first encoding stage consists of 

2 'P polar code encoders of length 2 ^, and the second encoding 
stage consists of 2'J polar code encoders of length 2^*. This is 
shown in Pig. |3] where we denote the polar encoding structures 
B 2 p{Ff^) and by and respectively. The 



Fig. 3. Encoding scheme based on ® 

PPA is run on the base code of length 2^ to find the puncturing 
sequence. To obtain a code rate of k/{N — m2'^), the first m 
bits in the progressive puncturing sequence are punctured at 
the output of each code at the second encoding stage. 
This results in a regular puncturing pattern on the long code, 
derived from base code’s progressive puncturing pattern. This 
regular puncturing pattern results in having the output channels 
of each G®'* code be transmitted on multiplicities of the same 
channel type, either punctured channels with zero capacity or 
the physical transmission channels. Hence, this construction 
is a generalization of the compound polar code construction 
0, and cotTesponds to an 2^ multi-channel compound polar 
code, which we will show is beneficial when mapping the 
bits to symbols of higher order modulations. To satisfy HARQ 
requirements, the information set is fixed across all codes with 
different rates in the same family with mother code length N, 
and is chosen to minimize the BLER of a high rate code 
in the RC family. Besides having the best performance at 
high code rates, selecting the information set at the high-rate 
code guarantees that no zero-capacity channels exist in the 
information set of the lower-rate codes. 

C. Rate-matching, channel interleaving, and bit-mapping 

In Section IIII-BI we proposed an approach to obtain a 
regular puncturing pattern of a large block length based 
on the puncturing pattern of smaller length selected by the 
progressive puncturing algorithm. However, with this regular 
puncturing structure, the punctured codeword length is re¬ 
stricted to £2"? and the transmission code rate is also restricted 
to k/{i2‘>), for £ G {0,1,..., 2^}. To achieve a finer resolution 
of code rates, we propose a simple rate-matching scheme that 
preserves the regular puncturing obtained by the two step 
polarization and the puncturing order on the base code found 
by the PPA. 

Eor a polar code of block length N = 2", we follow the 
approach in Section IIII-BI and consider the encoding structure 











































write codeword bits 
into matrix form 


with two stages based on ® where p and 

q are positive integers and p + q = n. First, the 2" coded 
bits, (xi,X 2 t ■ ■ ^xm), are arranged into a 2"? x 2^ matrix 
such that the (i, j)th element of this matrix is X(i_i) 2 P+j, for 
1 ^ ^ 2'J and 1 ^ j ^ 2^’. The first row of this matrix 

is (xi, X 2 ,..., X 2 p). In this structure, each row corresponds 
to the component code, and each column corresponds to 
the component G®'^ code. Next, the 2^ columns are permuted 
according to the reverse of the puncturing order found by the 
PPA on the base code of length 2^. To be more specific, we 
rearrange the columns such that the codeword bits located in 
the last m columns of the permuted matrix are the punctured 
bits of the regular puncturing pattern which punctures m bits 
in every 2^ bits, for m € {1,2,...,2?’}. This rate-matching 
structure inherently does row-column channel interleaving to 
the polar codeword, where the transmitted bits are written row¬ 
wise into the matrix and read column-wise after the described 
column permutation. Let the desired transmission rate be k/L. 
The transmitted codeword bits, denoted as (xi, X 2 , • ■ • > 
are read from the permuted matrix such that X(j_i)29+i is the 
(i, j)th element of the permuted matrix. If L = N — m2^, then 
the last m columns will not be transmitted, and the two-step 
compound polarization with regular puncturing is preserved. 
Moreover, this rate-matching structure does not constrain L 
to be a multiple of 2'?, where the first L bits are always 
read column-wise, after column interleaving. From the PPA, 
the column interleaving guarantees that the non-transmitted 
(punctured) bits will have the least impact on the BLER given 
by dTji. Moreover, if L > the matrix is treated as a circular 
array, where the (TV -f l)th transmitted bit is the (1, l)th bit 
of the array, while reading column-wise. According to the 
application, further polynomial interleaving can be applied on 
the output interleaved sequence. 

For a codeword of length 2^^, the rate-matching structure is 
depicted in Fig|4] In this example, we consider the encoding 
structure based on i? 2 ^^ (F®® (g)F®^). The codeword of length 
2^^, (xi, X 2 ,..., X 4096 ), is arranged into an 2^ x 2® array. 
Then, the columns of the array are permuted according to 
the reverse of the puncturing order found by the progressive 
puncturing algorithm. The transmitted codeword is read from 
this permuted array column-wise. 

In case bit-interleaved coded modulation (BICM), this 
structure can be combined with bit-mapping to higher-order 
modulated symbols, such as 16-QAM or 64-QAM. We apply 
a particular bit-mapping scheme under this rate-matching 
structure to satisfy the compound polarization structure. The 
transmission channel with M-QAM modulation and gray 
mapping can be modeled as a multi-channel with log 2 (M )/2 
different channel types ||7], where BICM detection results 
in log2(M)/2 BICM channels of different reliability. The 
proposed bit-mapping is based on our generalized compound 
polar code structure, where the 2 ^ multi-channels of the base 
polar code are mapped to the log2(M)/2 BICM channels and 
the zero-capacity channels resulting from puncturing. 

An example for the bit-mapping from the binary coded bits 
to the 16-QAM symbols is shown in Fig. |5] where the two 
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Fig. 4. Example of rearrangement of codeword of length 2^^ with the rate¬ 
matching structure 
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Fig. 5. Bit-mapping from codeword bits to 16-QAM symbols in matrix form 


constituent BICM channels from the 16-QAM modulation are 
denoted as Wi and W 2 - This is combined with rate-matching, 
where the the first half of transmitted bits are transmitted on 
the first constituent channel and the second half are transmitted 
on the other constituent channel. 

Punctured polar coding with 64-QAM modulation and gray 
mapping can be also modeled as polar coding over a multi¬ 
channel with four different channel types, three channels 
with different reliabilities are due to 64-QAM BICM with 
gray-mapping and the fourth channel is the zero-capacity 
channel due to puncturing. After rate-matching, the transmitted 
columns are grouped into three groups of the same size and 
each mapped to a different BICM channel type. The idea in 
this bit-to-symbol mapping is to guarantee, whenever possible, 
that each component G®"^ polarizes over the same channel 
type. Hence, different assignments of the columns to the BICM 
bit-channels are possible. 

D. HARQ Chase combining and incremental redundancy 

The proposed rate-matching scheme allows transmission of 
rate-compatible polar codes with both Chase combing (CC) 
and incremental redundancy (IR). If the decoding result of 
the previous transmission failed, then the same modulated 
symbols are re-transmitted with CC, but different symbols 
can be transmitted with IR. Assume the transmission rate of 
k/L is fixed across all HARQ transmissions and a maximum 











































of t transmissions. The proposed structure allows seamless 
IR or Chase transmissions at any rate. With HARQ CC, 
all transmissions will send the L bits column-wise from the 
permuted matrix, starting from the first column. With IR, 
the rth transmission may start sending the L bits column¬ 
wise starting from ((r — l)2P/t + l)th column. The BICM 
bit-mapping is then applied by circularly shifting the BICM 
column assignments (r — l)2P/t columns. Other combinations 
of the IR column assignments of each transmission with 
the BICM bit-channel assignments that preserve the two-step 
polarization structure are also possible. 


IV. Analysis of proposed RC polar codes 


First, we provide a theoretical analysis on the achievable 
rates, and then we provide numerical simulation results. 

Consider the modified polar transformation of Fig. |3] with 
base polarization block G®^, N' = 2^, and a puncturing 
pattern V of size m. Instead of being transmitted on iden¬ 
tical channels as in Arikan’s construction, the punctured bits, 
Xi : i € V, ctre transmitted on zero capacity channels and the 
remaining XiS are transmitted through N' — m independent 
copies of W. Let yi = 0 for i G V. The modified channel 
Wn', between the encoded information sequence and 
the channel output , are described by the transitional 
probabilities WN'{yi \u^') = where 

[A^'] is the set of indices 1,2,..., N'. Then, the modified bit- 
channels for i G [W'] are given by 




E 




( 8 ) 


The sum-capacity can be found by generalizing jT) lemma 4], 
where C{W) denotes the symmetric capacity of channel W. 


Lemma 1. For any m < N', N' = 2^, and any puncturing 
pattern V on the output of G®^ with cardinality m, the sum 
capacity is ~ Tn)C{W). 


Flence, the portion of almost noiseless channels approach 
the average of the capacities of which is equal to 

(1 —^)C(kF) by Lemma[T] Applying the channel polarization 
theorem 12 and Q Theorem 5], we show that the two-step 
polarization of Fig. |3] results in a family of capacity-achieving 
punctured polar codes. 


Theorem 1. Consider a polarization base of length N' = 2^ 
with m punctured bits and j3 < 0.5. There exists a family 
of two-step polarized codes approaching the rate Cm = 
(1 — ^)C{W) with vanishing probability of error bounded 
by 2~^^, as code length N goes to infinity. 


Note that Cm is the average channel capacity after punctur¬ 
ing a fraction of m/N' output bits. The above theorem can be 
generalized to Cm = (1 ~ ;^)E[C(kF)] for the case of BICM 
channels, where E[C(kF)] is the average symmetric capacity 
of the BICM channels as described in subsection Illl-CI 
Next, we provide numerical simulations for the proposed 
interleaved RC polar codes with BICM. The progressive 
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Fig. 6. BER of RC polai' codes with a fixed information set of size 88 at 
different code rates under AWGN (top) and fast fading (bottom) channels 


puncturing sequence is found on a base code of length 32, from 
which two RC families with mother code lengths N = 256 and 
N = 1024 are constructed. Let sj, be the fcth transmitted QAM 
symbol, then the fcth received symbol is = akSk+rik, where 
Ofe is the zero-mean complex Gaussian fading coefficient or the 
Rayleigh fading coefficient for QAM or BPSK, respectively. 
Uk is the additive white Gaussian noise (AWGN) with signal 
to noise ratio (SNR) 10log]^Q(l/(2cr^) in case of QAM, or 
101og]^Q(l/tT^) in case of BPSK. The puncturing sequence 
(0, 16, 8, 24, 2, 20, 26, 12, 10, 18, 4, 22, 25, 6, 13, 14, 
1, 17, 28, 3, 5, 9, 29, 11, 19, 7, 21, 15, 23, 27, 30, 31) 
is found by running the PPA with the GA at SNR = 3.5 dB 
on the base code with rate 11/32. 

Fig. |6] shows the bit error rate (BER) performance by SC 
decoding of the RC family with N = 256 and \I\ = 88 bits at 
code rates R G {0.9,0.8,0.7,0.6,0.5,0.4,11/32} for BPSK 
transmissions on both AWGN and fast fading channels. To 
select the information set, the bit-channel error probabilities 
are estimated by genie-aided SC decoding of the codes with 
rate R = 88/98 at SNRs of 3.5 dB and 6.5 dB on AWGN 
and fading channels, respectively. 

In Fig. |7] the normalized throughput with IR and CC HARQ 
schemes and BICM is illustrated for the RC family with 
N = 1024 and \X\ = 352 bits. A modulation and coding 



















scheme (MCS) is defined by the code of rate R and the 
modulation order M. The normalized throughput is defined 
as T = i?log 2 (M)(l — BLER)/f, where BLER is the average 
residual block error rate after all decoding attempts and t is 
the average number of transmissions required for successful 
decoding of an information block at the tested SNR. Assuming 
a maximum of f = 4 transmissions on an 16-QAM AWGN 
channel, with 3 RC polar codes with rates from 0.34 to 
0.92. It is observed that IR provides about 3 dB gain over 
CC. By comparing IR and CC at i? = 11/32, we conclude 
that changing the bit-to-symbol mappings across transmissions 
provides extra diversity and better performance. In system 
design, one selects the MCS with the highest throughput at 
a given SNR. To maximize throughput, it is best to utilize 
16-QAM with R = 11/32 till an SNR of 1.2 dB, then use 
a transmission rate of i? = 11/16 till an SNR of 12.8 dB, 
at which it is best to switch to 64-QAM with R = 11/16 
up to SNR of 19.7 dB, then use the highest transmission rate 
R = 11/12 with 64-QAM at higher SNRs. 

V. Conclusion 

To construct families of rate-compatible polar codes, a low- 
complexity algorithm that progressively selects the puncturing 
order on a base code of short length is devised and shown to 
have near-optimal performance. Based on a two-step polariza¬ 
tion construction, a flexible rate-matching scheme selects the 
transmitted bits for HARQ transmissions. Whereas previous 
algorithms would require storing the puncturing patterns for 
each designed code rate and length, the proposed scheme 
only needs to store one short (e.g. of length 32) puncturing 
sequence, which will be used to generate the puncturing 
pattern at any desired code rate and length. The proposed 
rate-matching and bit-mapping scheme preserves the code 
polarization in the presence of punctured bits and multi¬ 
channel bit-interleaved coded bits, and can result in any code 
rate by either puncturing or repetition of the coded bits in a 
particular order. IR and CC HARQ schemes with adaptive bit- 
interleaved polar-coded modulation on wireless channels were 
investigated. It is shown in this paper that punctured polar 
codes can have a capacity achieving property. 

Eurthermore, the proposed scheme can be applied with 
other polar code constructions such as the reduced-complexity 
relaxed polar code constructions ||9l or concatenated construc¬ 
tions that have improved burst-error correction capability ifTOll . 
With concatenated decoding IfTOll or list decoding IfTTI . the 
performance of RC polar codes is expected to be comparable 
to those of RC turbo codes and RC LDPC codes, cf. ifT^ . 
without suffering from error floors. The polar code rate¬ 
matching and bit-mapping schemes of this paper make the 
adoption of polar codes in future wireless systems more 
practical. 
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